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ABSTRACT 

We explore the stellar metallicity distribution function of the Galactic halo based on SDSS ugriz photometry. 
A set of stellar isochrones is calibrated using observations of several star clusters and validated by comparisons 
with medium-resolution spectroscopic values over a wide range of metal abundance. We estimate distances 
and metallicities for individual main-sequence stars in the multiply scanned SDSS Stripe 82, at heliocentric 
distances in the range 5-8 kpc and \b\ > 35°, and find that the in situ photo metric metallici ty distribution has 
a shape that matches that of the kinematically-selected local halo stars from lRvan & Norrisl We also examine 
independent kinematic information from proper-motion measurements for high Galactic latitude stars in our 
sample. We find that stars with retrograde rotation in the rest frame of the Galaxy are ge nerally more m etal 
poor than those exhibiting prograde rotation, which is consistent with earlier arguments by Carol lo et al.l that 
the halo system comprises at least two spatially overlapping components with differing metallicity, kine matics, 
and spatial distributions. The observed photometric metallicit y distribution and that of iRvan & Norrisl can be 
described by a simple chemical evolution model bv lHartwickl (or by a single Gaussian distribution); however, 
the suggestive metallicity-kinematic correlation contradicts the basic assumption in this model that the Milky 
Way halo consists primarily of a single stellar population. When the observed metallicity distribution is de- 
convolved using two Gaussian components with peaks at [Fe/H] pa —1.7 and -2.3, the metal-poor component 
accounts for ~ 20%-35% of the entire halo population in this distance range. 

Subject headings: Galaxy: abundances — Galaxy: evolution — Galaxy: formation — Galaxy: halo — Galaxy: 
stellar content 



1. INTRODUCTION 

Knowledge of the nature of the stellar Galactic halo, 
which collectively preserves a detailed rec ord of our Galaxy's 
formation in the e arly universe (e.g., Eg g en et"afl 119621: 
ISearle & Zinnl fl978). has expanded a great deal in the past 
few years. With the advent of large photomet ric surveys 
such as the Sloan Digital Sky Survey (SPSS: lYork et all 
2000; Stought onetalJBOOa lAbazaiian et ail 120031 12004 , 
20051 lAdelman-McCarthv et al I 120061; iGunn et al.l 2006 ; 
Adelman-McCarthv et alj|2007l 120081; lAbazaiian et al.ll2009 ; 
Aihar a et al. 1 120111) . as well as the massive spectroscopic 
follow-up efforts that have come from SDSS, in particular, 
the Sloan Extens ion for Galactic U nderstanding and Explo- 
ration (SEGUE; lYannv et"aT] 120091) . and the recently com- 
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pleted SEGUE-2 extension (Rockosi et al., in preparation), 
the opportunity to collect a vast amount of information on 
the nature of the Galactic halo has arrived. Detailed quantita- 
tive comparisons of these observations with the spatial, kine- 
matical, and chemical space inferred from theoretical predic- 
tions of models for the formation of the stellar halo through 
the hierarchical merging of higher- mass subhalos and accre- 
tion of lower-mass subhalos (e.g..lBullock & Johnstonl | 2005 : 
Johnston et all 120081; IPont et all 1201 It iMcCarthv etai1l2012 : 
Tissera et al. 2012) are now within reach. 

Much attention has been paid to understanding the nature 
of the stellar halo(s), in terms of the chemical and kinemat- 
ical properties of its constituent stars. For example, based 
on the medium-resolu tion SDSS spectroscopy of "local" halo 
stars (c/heiio < 4 kpc), iCarollo et all (120071 120101) argued that 
our Milky Way stellar halo is a superposition of two over- 
lapping systems, the inner and outer halos, that are distinct 
in metallicity, kin ematics, and spa tial distributions (see also 
iBeers et al.l 120121) . Furthermore, iNissen & Schuster! (120101) 
showed that high [a/Fe], metal -poor halo stars in the solar 
neighborhood are mainly on prograde orbits, while those with 
low [a/Fe] abundances are preferentially found on retrograde 
orbits. High [a/Fe] halo stars are likely associated with a dis- 
sipative component of the Galaxy that experienced a rapid 
chemical evolution, while the low [a/Fe] stars could be ac- 
creted components from dwarf galaxies that had lower star- 
formation rates. For example, one could point toward low- 
mass (low metallicity) dwarf-like galaxies (surviving exam- 
ples of which might include the ultra-faint dwarfs discovered 
by SDSS, e.g., Willman et al. 2005, Belokurov et al. 2006a,b, 
Zucker et al. 2006, and many others since; see also the dis- 
cussion of Carollo et al. 2012 and Qian & Wasserburg 2012) 
as likely progenitors of the outer-halo stellar population of the 
Milky Way. 
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The key observational ingredient in testing and studying 
the dual nature of the Galactic halo is the metal abundance 
for individual stars. In particular, such a diversity of stel- 
lar populations in the halo could have been imprinted on 
the metallicity distribution function (MDF) of halo stars. 
Previous results on the halo MDFs have been based pri- 
marily on spectroscopic surveys of likely halo stars se- 
lected from k inematically-biased (e.g jRvan & Norrisll991bt 
ICarnev et alj|1996l) or kinematica lly-unbiased [e.g., the HK 
survey of lBeers et all dl985l [1991 : t he Ham burg/ESO survey 
described by iChristliebl ( 120031) and iChristlieb etail (120081) 1 
searches for metal-poor stars. These approaches have the ad- 
vantage of reasonably efficiently identifying stars of the low- 
est metallicity, which are of great interest in their own right, 
but they suffer from selection biases that can be difficult to 
quantify. Biases are of particular importance if one seeks to 
understand the global properties of the stellar halo, a neces- 
sary step in order to tell a coherent story of the assembly and 
evolution of large spiral galaxies such as the Milky Way. In 
this regard, large photometric surveys such as SDSS can be 
used to recover relatively unbiased information on the nature 
of the Galactic halo. 

Broadband photometry can provide reasonably accurate es- 
timates of stellar metallicities, temperatures, and distances, 
at least for main-sequence stars. In the pioneering work 
to recover stellar m etallicities from SDSS ugriz photometry, 
Ivezic et al. (2008a) devised a method of using polynomial re- 
gressions based on spectrosco pic calibrations of de-reddened 
u-g and g— r colors (see also Peng et al.l l2012l) . which is sim- 
ilar to the traditional method of relating U V excess wit h posi- 
tions on the Johnson UBV diagrams (e.g.. lCarnevl[T979l) . The 
clear advantage of using a photometric metallicity technique 
is the efficiency of estimating metallicities for main-sequence 
stars, which are the most plentiful and representative sample 
of individual stellar p opulations. 

As recognized by llvezic etafl d2008al) . the polynomial- 
based photometric metallicity technique they developed be- 
comes quite insensitive at lower metallicity. They compared 
photometrically estimated metallicities with those obtained 
for a subset of the same sample of stars with available spectro- 
scopic determinations, and found that the photometric [Fe/H] 
estimates effectively saturate at [Fe/H] spec < -2, and cannot 
be extended to lower metallicities due to limits placed by the 
level of photometric errors obtained by SDSS. 

Although the power of photometric metallicity determina- 
tions from broadband colors is fundamentally limited at low 
metallicities, metallicity estimates can indeed be pushed down 
to below [Fe/H]^ -2.5, through the use of well-defined c olor- 
metallicity relations. For example, in I An et al.l d2009bl) . we 
used theoretical isochrones, calibrated using sets of observed 
cluster data, to estimate photometric metallicities for a large 
number of stars in the halo, based on SDSS gri photometry. 
Although individual metallicity estimates exhibit significant 
errors, due to the weak metallicity sensitivity of gri colors, 
it was shown that application of the calibrated isochrones can 
extend the low metallicity limit down to at least [Fe/H]~ -2.5, 
based on a comparison with SDSS spectroscopic abundance 
measurements. 

In this paper, we extend our previous efforts on the cali- 
bration of stellar isochrones to provide improved photomet- 
ric metallicity estimates, taking into account the full set of 
information from the SDSS ugriz measurements (§ |2). A 
newly calibrated set of isochrone s is dev eloped, supersed- 
ing the set discussed by lAn et al.l (l2009al) . Based on these 



new isochrones, we re-explore application of the photomet- 
ric metallicity estimation technique [ which is clearly indep en- 
dent of that discussed previously by llvezic et ail d2008al) l to 
samples of stars with available SDSS photometry (§|3}, in or- 
der to construct an unbiased MDF of the Galactic halo (§|Ui. 
Based on our derived photometric MDFs and a limited in- 
vestigation of available kinematic information, we show that 
simple chemical evolution models that are assumed to apply 
to a single stellar population provide inadequate descriptions 
of the nature of the Galactic halo. 

2. PHOTOMETRIC METALLICITY TECHNIQUE 

The key ingredient for photometric estimation of stel- 
lar metallicities (as well as distances) is the use of well- 
established stellar color-magnitude relations over a wide 
range of stellar atmosphe ric param e ters. Fo r this purpose, we 
adopted the f ormalism in | An et al.l d2009al) . where we used a 
set of YREC dSills et al. 2000) stellar isochrones in the SDSS 
filter system calibrated ag ainst observed u griz (open and glob- 
ular) cluster photometry dAn et al.l l2008) to estimate the stel- 
lar metallicities and distances of individual stars. 

As an initial exploration of this technique, we used the 
SDSS gri data to constrain (median) photometric metallici- 
ties ( hereafter [Fe/Hl phnt) for bulk stellar populations in the 
halo dAn et al.ll2009bl) . However, it was necessary to include 
M-band measurements to tightly constrain photospheric metal 
abundances, and therefore to construct a precise MDF. This 
has a consequence of limiting the application of the photo- 
metric technique to nearby main-sequence stars, because of 
the relatively shallow survey limit in the SDSS u passband 
(u = 22.0 at the 95% detection repeatability for point sources). 

In this section, we describe the development of empirical 
corrections to color- T e s rel ations for stellar isochrones in the 
native SDSS filter system dAn et al.ll2009al) . including the u 
passband, and test the accuracy of color- r e ft relations and pho- 
tometric metallicities using spectroscopic estimates from t he 
SEGUE Stellar Parameter Pipeline (SSPP:|Lee et al.ll2008allbfc 
lAllende Prieto et al.ll2008HSmolinski et al.ll201 flT~ 

2.1. Background: Calibration of Isochrones 

In the following analysis we employ the same set 
of underlying stellar interior mod els and t he sam e a- 
element enhancement scheme as in lAn et al.l (l2009al) . mo- 
tivated by the observed behavior of these elements among 
field dwarfs and cluster stars from spect roscopic studies 
(e.g., IVenn etalJ 120041 iKirby et al.l 12008): [a/Fe]= +0.4 
at [Fe/H]=-3.0, [a/Fe]= +0.3 at [Fe/H]= -2.0,-1.5,-1.0, 
[a/Fe]= +0.2 at [Fe/H]= -0.5, and [a/Fe]= +0.0 at [Fe/H]= 
-0.3,-0.2,-0.1,+0.0,+0.1,+0.2,+0.4. A linear interpolation 
was made in this metallicity grid to obtain isochrones at in- 
termediate [Fe/H] values. We adopted an age of 13 Gyr 
for -3.0 < [Fe/H] < -1.2, and 4 Gyr at -0.3 < [Fe/H] < 
+0.4, with a linear interpolation between these two metallicity 
ranges. As described below, our strategy minimizes the appli- 
cation of our calibrated stellar models to main-sequence turn- 
off stars, for which colors are most sensitive to the adopted 
age. For this reason, our main results are insensitive to the 
age assumption, but readers are cautioned if they apply our 
models to stellar populations with different ages. The above 
age-metallicity-[a/Fe] relations are adopted throughout this 
paper. 

Colors and magnitudes predicted by theoretical models typ- 
ically do not agree with those observed from star clusters. 
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Therefore, we define d color- r e ff- [Fe/H ] corrections to better 
match the data. In lAn et al.l d2009al) . we used M67 clus- 
ter data to define color- T e s relations in u — g, g — r, g — i, 
an d g — z, at the cl uster's metal abundance ([Fe/H] = 0.0). 
In lAn et al.l (1200 9b). we adopted these correction factors in 
gri color indices at solar abundance, and used a linear ramp 
between [Fe/H]= -0.8 and [Fe/H]= 0.0. In this weighting 
scheme, empirical corrections become zero at [Fe/H]= —0.8 
and below, so metallicity estimation of metal-poor stars are 
essentially those obtained from pure theoretical calculations, 
and are not affected by the empirical color corrections. This 
choice was motivated by the fact that the models are in satis- 
factorily agreement with data in gri color indices for globular 
clusters dAn et al . 2009a). Above solar abundance, we applied 
the M67-based color corrections to all models, assuming that 
offsets in colors between the model and data are independent 
of the metal abundance of a star. 

2.2. Updates: Calibration oflsochrones 

In principle, cluster sequences alone can be used to directly 
estimate stellar distances and metallicities, as long as there 
exist well-measured cluster-star samples over a wide range of 
[Fe/H]. However, cluster data are often noisy (e.g., some clus- 
ter fiducials are found to be essentially superimposed, even if 
their [Fe/H] values are different), so we choose instead to be 
guided by theory, in order to infer stellar colors at any given 
metallicity (and/or age). 

In particular, updates on the empirical calibrations of mod- 
els were necessary to adjust model u — g colors, which have 
a profound impact on the accuracy of photometric metallicity 
estimates, to better match observed cluster fiducial sequences 
over a wide range of metal abundances. Below we describe 
our adopted methodology used to generate a new set of cali- 
brated isochrones for all of the ugriz color ind i ces. 

We expanded upon the strategy of lAn et aTl (l2009al) for the 
color-r e ff calibration by t he addition of cluster fiducial se- 
quences in lAn et al.l (120081) over a wide range of metal abun- 
dances. Although our pure theoretical models predict colors 
that are consistent with observed fidu cial sequences within 
the total systematic and random errors (An et al. 2009a), there 
still remains a small, but suggestive systematic residual pat- 
tern of the color offsets over T e ff and [Fe/H]. In the updated 
calibration set, we have attempted to minimize these effects 
to better constrain estimates of the stellar parameters. 

In the color- r e ff- [Fe/H] calibration, which is described in 
detail below, we used cluster fid ucial sequences for sev- 
eral globular and open clusters in lAn et al.l (120081) . These 
clusters are listed in Table |2] along with our adopted val- 
ues for the cluster parameters, which are t he same as thos e 
used in our earlier model comparisons (lAn et al] [2008). 
The [Fe/ H] and E(B - V) es timates for globular clusters 
are from iRraft & Ivansl (120031) . who used Fe II lines from 
high-resolution spectra to compare colors derived from high- 
resolution spectroscopic determinations of T e s with the ob- 
served colors of the same stars. Distances to the glob- 
ular clusters are all Hi pparcos-based subdw arf fitting dis- 
tances. We adopted the ICarretta et al] (120001) distance esti- 
mates whenever they are available, since they employed the 
same metallicity scale for both subdwarfs and cluster stars 
in the su bdwarf-fitting t e chniq ue; otherwise, we adopted dis- 
tances in IRraft & Ivansl (120031) . For M67, we took the aver- 
age reddening and metallicity es timates from hig h-resolution 
spectroscopy in the literature dAn et al] 12007b , and refer- 
ences therein), and adopted a cluster distance estimated from 



TABLE 1 

Photometric Zero-Point Differences 



Cluster 
Name 


SDSS 
Runs 




UberCal - Photo 




{Au) 








<Az) 


M15 


ZD 66 


0.020 


0.007 


-0.005 


-0.006 


—0.010 


M15 


1739 


0.008 


0.003 


-0.005 


0.003 


—0.009 


M92 


4682 


-0.037 


-0.013 


-0.010 


-0.015 


-0.034 


M92 


5327 


-0.001 


-0.004 


-0.031 


-0.035 


-0.026 


M13 


3225 


0.011 


-0.003 


-0.001 


-0.005 


0.006 


M13 


3226 


0.007 


-0.011 


-0.006 


-0.002 


0.007 


M3 


4646 


0.039 


0.014 


0.016 


0.010 


0.008 


M3 


4649 


-0.022 


0.011 


0.013 


0.004 


-0.015 


M5 


1458 


-0.002 


-0.005 


-0.001 


0.004 


-0.001 


M5 


2327 


-0.010 


-0.001 


-0.007 


0.008 


0.019 


M67 


5935 


-0.002 


0.001 


-0.002 


0.001 


-0.004 


M67 


5972 


0.004 


-0.001 


-0.001 


0.007 


-0.009 


NGC 6791 


5416 


0.003 


0.003 


0.004 


0.004 


0.008 


NGC 6791 


5403 


0.001 


0.001 


0.001 


0.000 


0.000 



an empirically calibrated set of isochrone s in the Johnson- 
Cousins-2MASS system dAn et al.ll2007bb . For NGC 6791, 
we adopted the average [Fe/H] f rom high-resolut ion spectro- 
scopic studies (see references in lAn et al]|200 9a^. The clus- 
ter's reddening and distance estimates are based on the appli- 
cation of our calibrated isochrones in the Johnson-Cousins- 
2MASS system (An et al., in preparation). 

We employed UberCa l (Uber-calibration) magnitudes 
(Padmanabhan et al. 2008) for calibrating the cluster systems, 
instead of the "Photometric Telescope (PT)"-calibrated mag- 
nitudes (hereaft er Photo mag nitudes)!!] Our original clus- 
ter sequences in I An et aTl (120081) were on the P hoto system, 
where the standard SDSS photometric pipeline dLupton et al] 
2002) was used to define stellar colors and magnitudes in 
SDSS. However, iPadmanabhan et al] d2008l) later devised a 
method of improving a relative photometric calibration error 
by using repeat measurements in the overlapping fields of the 
survey. Since DR7, SDSS takes the UberCal magnitudes as 
the default magnitudes. 

We transformed the fiducial sequences in the Photo sys- 
tem onto the UberCal system, by applying zero-point dif- 
ferences between the two systems in the cluster flanking 
fields, where we derived pho tometric zero points for the clus- 
ter fiducials dAn et al.l 12001) . The differ ences between the 
two calibrations are not alarmingly large for the fiducial se- 
quences used in the current study. The mean differences are 
+0.003 ±0.014 mag, +0.003 ±0.015 mag, -0.001 ±0.010 
mag, -0.005 ±0.028 mag, and -0.003 ±0.013 mag in ugriz, 
respectively, in the sense of UberCal minu s Photo magni- 
tudes, for all the cluster fields considered in lAn et all (120081) . 
except from the few imaging stripes from which we could not 
retrieve UberCal magnitudes. The zero-point differences 
between UberCal and Photo systems for several cluster 
fields used in this work are listed in Table Q] 

Gray cross points in Figures [T]|4]represent the color differ- 
ences between cluster fiducial sequences and pure theoretical 
stellar models in g - r, g-i, g-z, and u-g, respectively. In 
each of the color indices, we have arranged the comparisons 
such that comparisons for metal-poor clusters are shown in 
the upper panels, and those for metal-rich clusters are shown 
in the lower panels. A gray strip represents a ±1<t range 
of a total systematic error in the comparison, including er- 
rors from the distance, reddening, age, and an assumed pho- 

13 See also http : / / www .sdss3.org/dr8/ algorithms / f luxcal . php 
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FIG. 1 . — Updated color- T e ff-[Fe/H] calibrations of YREC isochrones. Gray 
cross points are color differences in g - r between the YREC models and 
cluster fiducial sequences. A gray strip represents a ±lcr range of a total 
systematic error in the comparison. Individual cluster cases are displayed in 
increasing order of metallicity from the upper to lower panels. Solid lines 
represent our derived color- r e ff- [Fe/H] corrections, based on all of the cluster 
comparisons shown above, where M92 and M67 comparisons are used as a 
baseline (see text). Color differences for individual stars are shown for M67. 



TABLE 2 
Adopted Cluster Parameters 



Cluster 


(m-M)o 


E(B-V) 


[Fe/H] 


Age a 




Name 


(mag) 


(mag) 


(dex) 


(Gyr) 


References 


M15 


15.25 


0.10 


-2.42 


12.6 


1 


M92 


14.64 


0.02 


-2.38 


12.6 


1.2 


M13 


14.38 


0.02 


-1.60 


12.6 


1.2 


M3 


15.02 


0.01 


-1.50 


12.6 


1 


M5 


14.46 


0.03 


-1.26 


12.6 


1.2 


M67 


9.61 


0.04 


+0.00 


3.5 


3 


NGC6791 


13.02 


0.10 


+0.40 


10.0 
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REFERENCES. — References for (m- M)o, E(B-V) and high 
resolution spectroscopic [Fe/H] values: (1) Kraft & Ivans (2003); (2) 
Subdwarf-fitting distances from Carretta et al. ( 2000); (3) lAn et all 
(2007b), and references therein; (4) An et al. (2012, in preparation), 
and references therein. 
a Adopted ages of clusters. 



tometric zero-p oint error on the fiducial sequences (see also 
lAn et alj|2008l for more details) . 

Solid lines in Figures [T]-[4] show empirical color correc- 
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FIG. 2. — Same as in Figureff] but for g — i. 

tions derived using cluster fiducial sequences. We based 
our calibration using cluster fiducial sequences for both M92 
and M67. In other words, color-T^g relations were defined 
with the observed m ain sequences of M92 at [Fe/H]= -2.4 



dKraft & Iv ans 2003), and those at 



with respect to M67 (see lAn et al 



Fe/H]= 0.0 were defined 



2007b). The choice of 



M92 and M67 for the base case was due to the fact that 
these clusters are well-studied and have reliably determined 
distances, foreg round reddening, and metallicity estimates 
dAn et al.ll2009al see also Table [2]). They both also have wide 
r e ff coverage (see Figures QJBll. 

To be consistent with our earlier color cal ibration with M67 , 
we used the same color corrections as in lAn et alJ d2009a). 
These are polynomial fits to the points over 4000 < r en -(K) < 
6000, and are expressed as follows: 

A(g-r) = 7.610- 19.91061+ 17.2796» 2 -4.983(9 3 , (1) 

A(g- z) = 6.541- 16.3130+ 13.461f5 2 - 3.675f9 3 , (2) 

A(g- z ) = 3.755-7.5956i+4.535(9 2 -0.6726i 3 , (3) 

A(u - g) = -7 .923 + 24. 1 3 1 9 - 23 .7276» 2 + 7 .6276> 3 , (4) 

where = r e ff/5040 K. These corrections are shown as solid 
lines in Figures [TJS1 (second panels from the bottom), and 
are in the sense that the above values should be added to 
the model colors. For M92, we used the color differences 
between the model and the cluster fiducial sequences as the 
baseline of the model calibration at the cluster's metal abun- 
dance (see second panels from the top in FiguresQJQ- 
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A quadratic relation in [Fe/H] was then used to parameter- 
ize the metallicity dependence of the color-temperature cor- 
rections, assuming the following functional form: 

Ai(T eS , [Fe/H]) = 5,(r eff ) + 0[Fe/H] + e[Fe/H] 2 , (5) 

where A,- represents color corrections in each of g — r, g- i, 
g—Z, or u-g, while Si represents color- T e fs corrections in 
Equations (d}-© in the i th color index at solar metallicity. 
Here, Q and are coefficients to be derived from the fit, 
where we used fiducial sequences for M15, M13, M3, M5, 
and NGC 6791. In fact, the problem is reduced to a single 
parameter fit, since our models were defined to match the ob- 
served main sequences of M67 and M92 at their metal abun- 
dances. The solid curves in Figures [T]-|4] show the resulting 
color corrections on the color- T e ff- [Fe/H] space. We extrap- 
olated color corrections below [Fe/H]= -2.4 using the best- 
fitting parameters obtained from the above equation. 

The r e ff vs. [Fe/H] range covered by our calibrating sam- 
ple clusters is sparse and non-uniform. In particular, hot stars 
(?eff ^ 6000 K) are not covered at or near solar abundance, 
while cool, metal-poor stars (r e ff < 5000 K) are outside of 
our calibration range. This is mainly because the clusters 
are found at varying distances from the Sun, while all of 
the cluster images were taken in drift-scan or time-delay-and- 
integrate (TDI) mode, with the same effective exposure time 
of 54.1 seconds per band. As a result, turn-off stars in M67 
were too bright and lower main-sequence stars in M92 were 
too faint for SDSS to obtain reliable magnitudes. The cluster 



age is an additional factor that affects the non-uniform cover- 
age in the stellar parameter space. 

In order to obtain full coverage of the color corrections in 
the r e ff vs. [Fe/H] plane, we assumed that the color difference 
from the model for M67/M92 at the high/low T e ff end extends 
and remains constant to r e ff = 7000 K/4000 K, respectively. 
We note that this reasonable, but arbitrary, extrapolation of 
the color- r e ff relations has minimal impact on the results of 
our analysis, because our calibrating cluster sample is in fact 
representative of the majority of stars detected in SDSS (in 
terms of r et Y and [Fe/H], among other parameters). The pur- 
pose of the extrapolation is to obtain stable color- T e g- [Fe/H] 
relations going forward. Note that our model calibration is 
valid for main-sequence stars only. 

Our adopted age of 13 Gyr at -3.0 < [Fe/H] < -1.2 in 
the color calibration (§ 12. It is justified by our earlier result 
that the main-sequence turn-off ages of our calibrating clus- 
ters are approximately 13 Gyr, when pure theoretical YREC 
models are dire ctly used in the age estimation (see Table 9 in 
An et al. 2009a). To a first approximation, a systematic error 
in the adopted age in the calibration results in a scale error in 
our color- r e ff corrections. Therefore, special attention should 
be paid when applying the models to stellar populations with 
different age-metallicity relations. 

2.3. Stellar Parameter Search 

We applied calibrated stellar isochrones to the observed 
ugriz magnitudes, and searched for the best-fitting stellar pa- 
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rameters — T e ff (or stellar mass), [Fe/H], and an absolute 
magnitude (or distance) — by minimizing the \ 2 of the fit, 
defined as 

C^obs./ ^model./) 



X 



(6) 



for each star over -3 < [Fe/H] <+0.4. Here, X b s .i andX mo dei,! 
are the observed and model magnitudes, respectively, in the 
i th bandpass. Stellar mass, [Fe/H], and distances were set as 
free parameters (i.e., three parameters and five data points for 
each star). This is equivalent to fitting a model to the ob- 
served spectral energy distribution in the wavelength versus 
flux space. 

We adopted foreground dust estimates by Schlegel et al. 
(119981). with th eoretical extinction coefficients given by 
lAnetal.ld2009al) : 

/?a(YREC) = 



E(B - V) 



= [4.858,3.708,2.709,2.083,1.513] 

(7) 

where A = u,g,r,i,z, respectively. These extinction coeffi- 
cients were computed using t heoretical spectral energy dis- 
tributions and the standard Carde lli et al.l (Q~989) extinction 
cu rve. Our values are in good agreement with those provided 
in iGirardi et al.l (120041) . The above valu es lie in between the 
default extinction co efficients in SPSS (TSchlegel et al.ll 19981 
SFD98) and those in lSchlafly & Finkbeineri(l201 lL SF1D: 

fl A (SFD98) = [5.155,3 .793 ,2.751 ,2.086, 1 .479] , (8) 

fi A (SFll)= [4.239,3.303,2.285,1.698,1.263]. (9) 

In terms of the u - g colors, which have a significant impact 
on photometric metallicity estimates, our E(u — g) = A u -A g 
coefficient is 18% smaller than the SFD98 value, but 19% 
larger than the Schlafly & Finkbeiner (201 1) value. These dif- 
ferences are clearly important, and we discuss their effect on 
our derived photometric MDFs in § |4]below. 

2.4. Comparisons with the SSPP Spectroscopic Estimates 

Photometric temperatures and metallicities in our approach 
are derived simultaneously, based on the observed ugriz pho- 
tometry. Each of these parameters are primarily constrained 
by different portions of the stellar spectral energy distribution; 
u — g colors are mostly responsible for the photometric metal- 
licities, while griz colors are sensitive to temperatures. Never- 
theless, there is a moderate level of correlation between them, 
so checking a photometric temperature scale is an important 
step toward obtaining accurate photometric metal abundances 
of stars. 

In IPinsonneault et al] (120121) . we verified the accuracy of 
our photometric T e ff estimation procedure using the most re- 
ce nt temperature scale from the Infrared Flux Method (IRFM) 
by Casagr ande et al.l (120101) . Photometric temperatures were 
estimated from gri z photometry in the Kepler Input Catalog 
dBrown et al.ll201lh using a calibrated set of isochrones - the 
same set of models used in the current work - then they were 
compared to the IRFM temperatures derived from 2MASS 
J - K s colors. Overall, we found good agreement between 
these two fundamental temperature scales ((Ar e ff) < ±40 K) 
at 4400 K < r RfF < 6000 K, where the IRFM scale is well- 
defined; see IPinsonneau lt et al] d2012l) for detailed discus- 
sions on this comparison. 

Here we compare our photometric temperature estimates 
with the spectroscopic estimates obtained from medium- 
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FIG. 5. — Comparison between photometric and spectroscopic (SSPP) tem- 
perature estimates for the SEGUE sample. The photometric estimates are 
based on all of the ugriz photometry. 



resolution SDSS/SEGUE spectra (Yannv et al .1120091) ; com- 
parisons are shown in Figure [5] In this comparison we in- 
cluded 5 1 , 999 stars from the intial sample of 162, 645 objects 
that satisfy the following selection criteria: 

• Sources are detected in all five bandpasses 

• xL-,„/v < 3.0, where yL- is a minimum y 2 defined in 
Equation (|6]l and v (= 2) is a degree of freedom 

• a,, < 0.03 mag 

• k>gg(YREC) > 4.15 

• S/N (SSPP) >20/1 

• logg ( sspp) > 3.5, 

where the sbuscript YREC and SSPP indicate parameters es- 
timated from the model isochrones and SSPP, respectively. 

The first two criteria above require a good fit to the data, 
and the limit on the M-band measurement errors exclude data 
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FIG. 6. — Comparison between photometric and spectroscopic (SSPP) 
metallicity estimates for the SEGUE sample. The photometric estimates are 
based on all of the ugriz photometry. Three curves in the top panel show a 
median and an interquartile range, respectively. 



points with less well-defined photometric metallicity esti- 
mates. The logg limit restricts the analysis to main-sequence 
dwarfs. The last two criteria select dwarfs with reliable spec- 
troscopic abundance measurements. We applied further cuts, 
based on colors (g — r > 0.2) and magnitudes (r > 14 mag), in 
order to remove stars with extremely blue colors or saturated 
brightness measurements. Most of the stars in the original 
SSPP sample were rejected due to their large w-band errors. 
Even if a more conservative cut on Xmin/^ < 10 were used, 
there would still remain a total of 58, 335 stars in the compar- 
ison. 

As shown in Figure|5] there is a constant T e g offset between 
the photometric and SSPP temperature scales over 5000 K 
< ?eff < 6500 K, where the maximum deviation is less than 
100 K in this temperature range. However, the good agree- 
ment between our temp erature scale and the IRFM scale 
dPinsonneault et all 12012) ensures that our temperature esti- 
mates are more or less closer to the fundamental temperature 



scaled The random scatter seen in Figure [5] is ^ 100 K for 
individual T e & estimates. 

Figure [6] shows a comparison between the spectroscopic 
(SSPP) metallicities from SDSS/SEGUE and our photometric 
estimates for the same set of stars, as in Figure|5] The central 
line in the top panel shows a median photometric metallicity 
for each 0.2 dex bin in the spectroscopic metallicity, and the 
other two lines represent the first and the third quartiles. There 
appears an extended metal-poor tail at a given SSPP [Fe/H] in 
Figure [6] in particular at low [Fe/H] values. This is because 
the metallicity sensitivity of stellar colors degrades at lower 
metal abundances, due to the non-zero photometric errors. In 
the next section, we model these profiles based on simulations 
of artificial stars. As will be shown, artificial star tests indicate 
that the median photometric metallicities (central line in Fig- 
ure|6]l are insensitive to photometric errors and the lower limit 
in the model grid at [Fe/H]= -3, and exhibit a maximum devi- 
ation of only A[Fe /H] ~ 0. 1 dex from the true [Fe/H] down to 
[Fe/H]= -2.8. The overall good agreement between our pho- 
tometric metallicities and the SSPP values shown in Figure|6] 
with a maximum deviation of less than 0.2 dex, confirms that 
our photometric metallicity scale is not far from the spectro- 
scopic scale down to [Fe/H]~ -2.5, perhaps as low as ~ -3. 
However, in order to perform a more stringent test of the ac- 
curacy of our photometric metallicity estimates, we require 
more comparison star samples below [Fe/H]^ -2.5. 

Because we included w-band measurements, which are 
more sensitive to metal abundances than the other SDSS 
bandpasses, the observed dispe rsion in F igure [6] is smaller 
than the equivalent Figure 2 in lAn et al.l (l2009bl) . which is 
based only upon gri. The la random scatter (estimated from 
the interquartile range divided by 1.349, assuming a normal 
distribution) is about 0.2-0.4 dex for individual stars. The 
[Fe/H] estimates in SDSS/SE GUE ar e precise to ~ 0.2 dex 
for individual stars dLee et al. 2008a b; lAllende Prieto et al. 
l2008t fSmolinsk i et al.l201 ll) . so this comparison indicates that 
the photometric metallicities can be as precise as ~ 0.3 dex 
per star, when good photometry is available. The down-turn 
of the mean trend at [Fe/H] ~ -0.8 in Figure|6]is likely caused 
by our ad hoc assumptions on the [a/Fe] and age as a function 
of [Fe/H] in the model. 

An analysis of likely member stars of Galactic open 
and globular clusters and high-resolution spectra of 
SDSS/SEGUE stars indicates that the SSPP exhibits a 
tendency to estimate [Fe/H] higher by about 0.25 dex for 
stars with [Fe/H] < -3.0, in particular for cool giants. 
Although the SSPP underestimates metallicities for stars 
with super solar abundances, the effect is less than 0.1 dex. 
Nevertheless, it should be kept in mind that the reliability of 
measuring photometric metallicities eventually requires a test 
again st high-r e solution spectroscopy. 

In lAn et al.l d2009bl) . we tested the accuracy of photomet- 
ric metallicities from the application of gri photometry, and 
found that the photometric technique systematically underes- 
timates [Fe/H] p hot in the low metallicity range (A[Fe/H] ~ 
0.3 dex at [Fe/H]~ -1.6. We speculated that the mismatch 
between the SSPP and photometric metallicities were likely 
caused by either unresolved binaries in the sample or a small 
zero-point offset in the model colors. Although the effect 
of unresolved binaries in the ph otom etric metallicities can- 
not be completely ignored (see § 12.51 below), small system- 
atic color differences can actually induce a systematic trend 

14 Further note that T c ff is a defined quantity in L = 47rK 2 oT c ff. 
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FIG. 7. — Comparisons between photometric and spectroscopic (SSPP) 
temperature (top panel) and metallicity estimates (bottom panel) for the 
SEGUE sample. Photometric estimates are based only on ugr photometry. 
Three curves in the bottom panel show a median and an interquartile range, 
respectively. 

in [Fe/H] p hot, especially when photometric metallicities are 
derived based on color indices that are weakly dep endent on 
Fe/H], such as the application of gri photometry in I An et ail 
2009b). The accuracy of the photometric technique (Fig- 
ure [6]l has improved as more metallicity-sensitive color in- 
dices (ugriz) are employed in the calculation of [Fe/H] p h ot . 

Figure|7]is the same as Figures [5}|6] but based only on ugr 
colors. The scatter in the comparison is larger than those seen 
in Figures |5]-|6] because of the weaker constraints on these 
parameters. The down-turn of the mean trend at [Fe/H] ~ -0.8 
is stronger than seen in Figure [6] Clearly, the inclusion of all 
ugriz passbands is preferred, in order to strongly constrain 
both photometric temperatures and metal abundances. 

2.5. Artificial Star Tests 



We performed artificial star tests in order to evaluate the 
effects of photometric errors and unresolved binaries and/or 
blends on our determinations of photometric [Fe/H] estimates. 
In particular, we use these simulation results to construct 
metallicity kernels, which are employed to assess the broad- 
ening of the [Fe/H] profiles resulting from the above effects. 
These kernels are then used in §|4]to deconvolve the observed 
MDF of the halo. 

Figure|8]shows the artificial star test results for a single stel- 
lar population without binar ies. In this exercise, we used the 
same set of models as in § 12.21 which include the main se- 
quence and a portion of the subgiant branch. For each of the 
various metallicity bins, we generated 20,000 artificial stars 
with masses above 0.65M Q from a model isochrone. The 
magnitudes were then convolved with a Gaussian having a 
standard deviation of 0.02 mag in gri, and 0.03 mag in the 
u and z passbands. We derived phot ome tric metallicities and 
distance estimates as described in § 12.31 the resulting MDFs 
are shown as solid black histograms in the left panels, and the 
distributions in distance modulus are in the right panels. We 
used the following selection criteria: 

• Xmin/^ < 3.0, where Xmin * s a minimum x 2 defined in 
Equation (O and v (= 2) is a degree of freedom 

• logg(YREC) > 4.15 

• 0.65 <M*/M sun < 0.75. 

The last criterion is specified in order to be consistent with the 
analysis of the actual data sets described in § |U and assumes 
that stars are on the main sequence. 

The solid red histograms shown in Figure [8] are those re- 
sulting from the simulations, but assuming photometric errors 
twice the size as above (0.04 mag error in gri, 0.06 mag error 
in // and z, and assuming no correlation between the errors in 
different bandpasses), illustrating the effect of the size of pho- 
tometric errors on the photometric metallicity estimates. For 
a given photometric error, the resulting dispersion is higher 
at lower metallicity, due to the weaker dependence of broad- 
band photometric co lors on metallicity at lower abundances; 
e.g., see Figure 1 in An et alj (I2009bl) . At solar abundance, 
the standard deviation of the [Fe/H] distribution in Figure|8]is 
0.1 dex, but increases to 0.3 dex at [Fe/H]=-1.6, and 0.5 dex 
at [Fe/H]= -2.4. At the lowest [Fe/H], the distribution be- 
comes asymmetric, with an extended low-metallicity tail. On 
the other hand, the distribution in distance modulus becomes 
progressively more symmetric at the lower input metallicity, 
because of the weaker metallicity sensitivity of colors and 
magnitudes at lower metal abundances. Our model set does 
not extend below [Fe/H]= -3 or above [Fe/H]= +0.4, which 
results in a piling up of stars at these metallicities. 

Figure [9] shows the effects of unresolved binaries and/or 
blends for three metallicity bins, [Fe/H]j nput = —1.6 (left pan- 
els), -2.0 (middle panels), and -2.4 (right panels), respec- 
tively. Photometric errors of 0.02 mag in gri and 0.03 mag 
in u and z passbands were assumed. To reduce the effects 
of interpolation errors in the isochrones, which sometimes 
produce a non-continuous distribution in [Fe/H] p h t, we con- 
volved the derived photometric metallicities with Gaussians 
having the grid size of the model (cr^e/H] = 0.1 dex). The 
top panels apply to stellar populations comprised entirely of 
single stars. The middle two panels show the results when un- 
resolved binaries are included in the sample. In total, 10,000 
primary stars were generated from a single [Fe/H] population, 
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FIG. 8. — Artificial star tests for various metallicity bins. Left: The thick black histogram shows a photometric MDF, where photometric errors of 0.02 mag in 
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[Fe/H]j nput values are shown on each panel, and indicated by the dotted vertical lines. No binaries are included in the simulation. Right: Same as in the left 
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and the same number of secondary stars was generated, based 
eithe r on a flat mass function, or adop ting the M35 mass func- 
tion dBarrado v Navascues et al.ll2001l) . These simulated stars 
were merged with 10,000 single stars with no secondary com- 
ponents to simulate a 50% unresolved binary fraction. The 
resulting magnitudes in ugriz were convolved with Gaussians 
having the specified photometric errors. 

As can be appreciated by inspection of Figure [9] the pres- 
ence of unresolved binaries has little impact on the derived 
photometric [Fe/H] distribution, because the dispersion is 
dominated by the effects of photometric errors. Binaries, due 
to the modified flux distributions due to a secondary com- 
ponent, typically have lower photometric metal abundances 
than those of the primaries alone, but these effects are mostly 
buried in the extended [Fe/H] distribution produced by the 
presence of non-zero photometric errors. The bottom pan- 
els in Figure [9] show photometric metallicity distributions of 
the SDSS/SEGUE spectroscopic sample in Figures [5}|6] in 



the same metallicity bins: -1.7 < [Fe/H]sspp < —1-5 (bot- 
tom left), -2.1 < [Fe/H]sspp < -1-9 (bottom middle), -2.5 < 
[Fe/H]sspp < -2.3 (bottom right). Except for a small offset 
in the metallicity scale between the SSPP and photometric 
metallicities (median [Fe/H] p h t values are sh own in the pan- 
els), which has already been pointed out in § 12.41 the overall 
[Fe/H] p hot distributions appear quite similar to those from the 
artificial star tests. 

Our strategy in this study is to deconvolve the observed 
MDF of halo stars using simulated profiles, and infer the un- 
derlying [Fe/H] distribution (§ HI). In order to carry out the 
deconvolution in a straightforward manner, we parameterized 
simulated [Fe/H] kernels by simultaneously fitting two Gaus- 
sian functions, as shown by the two green curves in each panel 
of Figure|9] The red curve is a sum of these two distributions. 
As discussed above, the precise binary fractions and the exact 
form of the adopted mass functions are not dominant factors 
in determining the overall shape of a given MDF. Therefore, 
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FIG. 9. — Effects of unresolved binaries and/or blends from artificial star tests at three input metallicity bins: [Fe/H]j npllt = —1.6 (left), —2.0 (middle), and -2.4 
(right). Top: A single star population without binaries. Middle: A 50% unresolved binary fraction with a flat, and a M35-type mass function, respectively, for 
secondaries. The green curves in each panel show Gaussian fits to the distribution; the red curve is their sum. Bottom: Photometric metallicity distributions of 
SDSS/SEGUE spectroscopic samples in the same metallicity ranges: —1.7 < [Fe/H]sspp <— 1.5 (left), -2.1 < [Fe/H]sspp <— 1.9 (middle), —2.5 < [Fe/H]sspp < 
-2.3 (right). Median [Fe/H] p i,ot values are shown in the panels. 



we opted to choose the M35 mass function with a 50% binary 
and/or blending fraction when comparing the observed MDF 
of the halo with simulated results. 

The top panel in Figure [10] displays photometric metallic- 
ity kernels for various input metallicity values ([Fe/H]; nput ), 
normalized to the peak value of each profile. These simu- 
lated profiles were generated from the same set of artificial 
star tests as above, assuming a M35-type mass function and 
by taking 0.02 mag error in gri and 0.03 mag error in u and z- 
The input [Fe/H] j nput values are from -1 .2 to -2.8 in intervals 
of 0.1 dex; for clarity, the simulated profiles are alternatively 
shown in thick and thin curves. Again, note that the resulting 
photometric distribution becomes broader at lower metallici- 
ties. 

The bottom panel in Figure [10] shows the median differ- 
ences between photometric metallicities and input [Fe/H] val- 
ues for the same data set used to construct the kernels in the 
top panel. Longer error bars represent the interquartile ranges, 
while the shorter error bars represent ±lcr errors, where we 
computed the ler error by dividing the interquartile ranges by 
1.349, assuming a normal distribution. 

3. SAMPLE SELECTION AND PROPERTIES 

Using the new photometric metallicity estimation technique 
described in the previous section, we estimated distances and 



metal abundances for individual stars in SDSS Stripe 82, 
which is one of the imaging stripes in SDSS that has been 
repeatedly scanned along the celestial equator. There are two 
photometric catalogs available in Stripe 82: the calibration 
or Standard Star catalog dlvezic et alj|2007l here after calibra- 
tion c atalogf^! and the coadd imaging catalog ( Annis et al. 
1201 ll hereafter coadded catalog). The calibration star cata- 
log contains stellar magnitudes for approximately one million 
sources, where the magnitudes were averaged at the catalog 
level. The coadded catalog is based on the coadded image 
products, and is about 0.5 mag deeper than the calibration cat- 
alog. Both catalogs, which were in principle produced from 
the same observations, obtained with the ARC 2.5-m SDSS 
survey telescope facilities, provide the most precise (~ 1%) 
photometry set available within SDSS, and therefore can be 
used to set the best available constraints on the photometric 
MDF of the Galaxy. Both catalogs were constructed on the 
Photo magnitude system. However, we directly employed 
our UberCal-based models in the parameter estimation us- 
ing this photometry, because the global photometric zero- 
point differences between the two systems are negligible. 

Below, we describe the selection of a p hotom etric sample 
from Stripe 82 designed to minimize bias (§ 13. U . and evaluate 

15 Available at http : //www . astro . Washington . edu/users/ivezic/sdss/ 
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FIG. 10. — Top: Normalized distributions of photometric metallicities at 
various input [Fe/H]j nput values. These simulated profiles were generated 
from artificial star tests at a given input [Fe/H]j npllt , assuming a 50% binary 
fraction and a M35-type mass function, and by taking 0.02-0.03 mag errors 
in ugriz (see text). Input [Fe/H] values are from —2.8 to —1.2 in an 0.1 dex 
increment, where simulated profiles are shown alternatively in thick and thin 
curves. Bottom: Median differences between photometric metallicities and 
input values. Longer error bars show interquartile ranges, and the shorter 
error bars represent ±l<x errors. The latter values were computed from the 
interquartile ranges divided by 1.349. 

the effects of unrecognized giants in the sample (§ 13.21) . In 
the subsequent section (§ |4), we present an unbiased MDF 
of the Galactic halo, and test the hypothesis that the halo is 
composed of two overlapping sub-components, based on both 
metallicity and kinematic information for our sample stars. 

3.1. Sample Selection and Bias 

An unbiased sample of stars is, of course, an important 
ingredient for obtaining a representative MDF of the Milky 
Way's halo population(s). Although photometric samples are 
less susceptible to sample biases than spectroscopic studies 
that make use of metallicity or color in their sample selection, 
there still exists a bias that needs to be taken into account, as 
discussed below. 

In this work we adopted a sample selection based on stellar 
mass, as estimated using our isochrones. FigurefTTIshows our 
color-calibrated models at [Fe/H]=— 2.4, — 1.6, — 1.2, and -0.8; 
the thick solid lines indicate where 0.65 < M*/M < 0.75. 
This range of stellar mass is similar to what is adopted in this 
work (see below). Our choice for the mass-based sample se- 
lection was motivated by our reasoning that, in order to obtain 
a representative sample of the halo, stars at different metallic- 
ities should be sampled in identical mass ranges. We consider 
this to be a superior choice to the more commonly adopted 
color-based selections, because of the strong relationship be- 
tween color, metallicity, and mass; narrow color-cuts in a stel- 
lar sample would produce a mix of stars, including less mas- 
sive, lower metallicity stars and more massive, higher metal- 
licity stars. Although colors can be used as a surrogate for 
temperatures, they have only a limited applicability for stel- 
lar masses, hence we believe our mass-based selection should 
produce a less biased sample of stars for the assembly of a 
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FIG. 1 1 . — The range of stellar mass in the model considered in this work. 
Color-calibrated models are shown at [Fe/H]= -2.4 (13 Gyr; left-most), -1.6 
(13 Gyr), —1.2 (13 Gyr), and -0.8 (9 Gyr; right-most), where soild lines 
represent 0.65 < M»/Mq < 0.75. The fiducial sequence of M13 is shown as 
a gray line for comparison. 



valid halo MDF. 

Figure [12] illustrates metallicity-luminosity relations as a 
function of stellar mass, where the luminosity is expressed in 
terms of a maximum heliocentric distance that can be reached 
by a star at a specific magnitude limit. Note that we adopted 
an age of 13 Gyr for models at -3.0 < [Fe/H] < —1.2, and 
4 Gyr at -0.3 < [Fe/H] < +0.4, with a linear interpolation 
between these two metallicity ranges (§ 12.11 ). The M-band 
magnitude limit was used in Figure [12J because of the strong 
sensitivity of this band on metal abundances. The w max = 20.6 
(top panel) and M max = 21.0 (bottom panel) correspond to a 
median photometry error of a « 0.03 mag in the calibration 
and the coadded catalogs, respectively. The err or si ze is simi- 
lar to those adopted in the artificial star tests (§ 12.51 ). 

At a given M-band magnitude limit, we computed a maxi- 
mum distance to which each star can be observed at a given 
stellar mass and metal abundance. At a fixed mass, metal- 
poor stars are brighter than metal-rich stars, thus they can be 
observed at greater distances than metal-rich stars (see Fig- 
urefTTTi. Note that a color-based selection would have the op- 
posite consequence on the luminosity of stars that would be 
included; at a fixed color (or temperature) stars are brighter at 
higher metallicity. It follows that, in a magnitude-limited sur- 
vey such as SDSS, color-based selection would produce sam- 
ples that are biased against more metal-poor stars at greater 
distances from the Sun. 

The gray shaded region in Figure [12] indicates the mass- 
distance limit set in our halo sample. In both panels, these 
areas are surrounded by a [Fe/H]= —1.2 model to insure that 
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FIG. 12. — Maximum heliocentric distance as a function of stellar mass that 
can be reached by a star at u < 20.6 mag (top panel) and u < 21.0 mag (bot- 
tom panel). These magnitude limits correspond to <r„ ss 0.03 mag error in u 
in the Stripe 82 calibration (top) and coadded (bottom) catalogs, respectively. 
The solid and dotted curves show distance limits at a number of different 
metallicity bins. The gray shaded region represents a mass-distance limit set 
for the halo sample in this work, which ensures that the sample is unbiased at 
[Fe /H] < — 1 .2, has a stellar mass less than the turn-off mass of the [Fe/H]= -3 
model, and is relatively free from thick-disk contamination (rfhelio > 5 kpc). 
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FIG. 13. — Distribution of Galactocentric distances for the final halo sam- 
ples in the calibration and the coadded catalogs in Stripe 82. 



the sample is unbiased at [Fe/H] < -1.2. Our photometric 
halo MDFs should be less affected by this choice, since we 
generally expect to obtain relatively few halo stars outside this 
metallicity range. An additional constraint on the stellar mass 
has been applied to the sample, because main-sequence turn- 
off masses (the high-mass end points of each curve in Fig- 
ure [12]) are varying at different metallicities. Therefore, we 



imposed an upper mass limit, M* < 0.75M Q , which is close 
to the main-sequence turn-off mass from a [Fe/H]= -3 model; 
turn-off masses occur at higher masses for [Fe/H]> -3. The 
lower limit to the heliocentric distance (5 kpc) was set to ex- 
clude possible thick-disk interlopers in the sample; this natu- 
rally results in a lower mass cut-off at ~ 0.65 M Q . We further 
impose d a c ut based on Galactic latitude (see discussion be- 
low in § 14. U . The application of these constraints results in a 
relatively narrow parameter space (gray areas in Figure [T2l> in 
the stellar mass versus distance plane (5 kpc < t/heiio < 8 kpc 
for the calibration catalog; 5 kpc < t/heiio ^ 9 kpc for the 
coadded catalog). In short, our sample selection (delineated 
by the gray regions in Figure [T2i i collects all of the stars at 
[Fe/H]< -1 .2 in a limited volume of the halo, given the mag- 
nitude limit set in a photometric catalog (cr„ < 0.03 mag in 
this work). The resulting distribution of our sample in Galac- 
tocentric distances is shown in Figure [T3] 

We have carried out comparisons of our samples obtained 
by both the calibration and the coadde d catalogs, using a 
search radius of 1". As recommended in llvezic et alj (2007), 
we only included sources in the calibration catalog with at 
least four repeated observations in each bandpass, for both 
the photometric comparisons and the following analysis, and 
proceed with the mean magnitudes and their standard errors. 
After carrying out iterative 3cr rejections, we found zero- 
point differences of -0.007 ±0.022 mag, -0.001 ±0.011 mag, 
+0.002 ±0.013 mag, and +0.003 ±0.015 mag, in u - g, g ■ - r, 
g - i, and g — z, respectively, for stars with u < 20.6 and fainter 
than 16 th magnitude i n each filter bandp ass; see also photo- 
metric comparisons in lAnnis et alj (1201 11) . The sense of these 
differences is that the u - g color measurements in the coad- 
ded catalog are redder than those in the calibration catalog. 
We restricted our comparison to only include those stars with 
Galactic latitudes at \b\ > 35°, to be consistent with the analy- 
sis of the halo MDF carried out below. The listed errors indi- 
cate the derived dispersions in color difference (rather than the 
errors in the mean), using ~ 100, 000 objects in each compar- 
ison. Because fainter M-band measurements yield less strong 
UV excesses, the use of photometry from the coadded catalog 
would lead to systematically higher photometric metallicity 
estimates than in the calibration catalog (this is confirmed in 
the next section). The rms deviations of the color differences 
between the two catalogs are stable (<7 co i or < 0.003 mag) along 
the 1 10° length of Stripe 82. 

Although a detailed study is beyond the scope of the present 
work, we identified a strong systematic deviation in the u 
passband (Am ~ 0.2 mag at u ~ 22 mag) at the faint end, be- 
yond the magnitude limit set in our sample (u < 20.6). One 
likely cause of the systematic offset is a Malmquist-type bias 
in the calibration catalog, because faint sources near the de- 
tection limit can either be detected or missed, due to the ex- 
istence of large Poisson errors. Since the calibration catalog 
took the average of the source magnitudes from individual im- 
ages, the mean magnitude could therefore be biased towards 
brighter source measurements. We avoided the photometry 
zero-point issue by requiring cr u < 0.03 mag; small photo- 
metric errors also make the photometric metallicity estimates 
more reliable. 

We have not made use of the main SDSS photometric 
database in this work. The 95% completeness limit in the 
u bandpass of the main survey is 22.0 mag, but the a„ = 0.03 
mag limit corresponds to u « 18.7 mag, on the order of 2 mag 
shallower than for the Stripe 82 catalog. At this limiting mag- 
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nitude, use of the main SDSS database would allow explo- 
ration of heliocentric distances only up to ^2.2 kpc for the 
construction of an unbiased [Fe/H] sample at [Fe/H]< —1.2, 
which is an insufficient volume to probe the halo MDF 

To summarize, in the following analysis we used the selec- 
tion criteria listed below: 

• Sources are detected in all five bandpasses 

• xliin/ 1 ' < 3.0, where x^in i s a minimum x 2 defined in 
Equation (O and v (= 2) is a degree of freedom 

• logg(YREC) >4.15 

• 0.65M Q < M* < 0.75M Q 

• Mass-metallicity-luminosity limits (gray region in Fig- 
ure [T2]i 

• <4eiio > 5 kpc 

• \b\ > 35°. 

The above selection restricts the photometric sample to 
0.2 < g - r < 0.5 (see FigurefTTli. The total numbers of objects 
that passed the above selection criteria are 2, 523 from the cal- 
ibration catalog and 2, 626 from the coadded catalog. The me- 
dian r magnitudes are 18.7 mag and 19.1 mag for the calibra- 
tion and the coadded catalogs, respectively, with a ~ 0.3 mag 
dispersion. The standard SDSS star-galaxy separation based 
on the difference between PSF (Point S pread Function) and 
model magnitudes is robust to r ~ 21.5 dLupton et al.ll2002t 
IScranton et al.l [2002), so the contaminati on by galaxies in 
our sample should be negligible (see also lAnnis et al.ll201 It 
iBovv et alj|2012ah . 

3.2. Contamination by Distant Giants and Thick Disk Stars 

Our calibration is valid for main-sequence stars only, and 
we have explicitly assumed that the great majority of stars 
in the sample are in their main-sequence phase of evolution. 
However, there certainly exist distant background giants and 
subgiants along each line of sight, and their distances can be 
greatly underestimated if a photometric parallax relation for 
main-sequence stars is blindly applied. Unfortunately, ugriz 
photometry alone is not sufficient to reliably disc riminate gi- 
ants f rom dwarfs, unlike horizontal-branch stars (llvezic et al.l 
120071) . and therefore some level of contamination by giants in 
our sample is unavoidable. 

We evaluated the level of expected giant and subgiant con- 
tamination using cluster photometry for both Ml 3 (rFe/H]= 
-1.6; see Figure [TT) and M92 ([Fe/H] = -2.4) (lAn et al.1 
2008). In order to carry this out, we applie d cuts based the 
o n y and sharp p arameters in DAOPHOT ( Stetsonl [l987l) as 
in lAn et all (12008b . Based on the large extent of these clusters 
on the SDSS CCD chips, we opted not to apply cuts based on 
the distance from the cluster centers. 

We assumed that the number density of the halo follows a 
power-law profi le with an index of -2.8, and a halo elliptic - 
ity of 0.6 (e.g.. lJuric et al.l 120081) . Using this model, along 
with the photometry for each cluster, we simulated color- 
magnitude diagrams for each line of sight in Stripe 82, and 
applied the same mass-metallicity-luminosity cuts to the sam- 
ple as in the previous section. We integrated the number of 
giants out to 30 kpc from the Sun, beyond which the giant 
contamination is negligible, because either those distant gi- 
ants are too faint to be included in our sample or they do not 



satisfy our selection criteria. Each object was tagged as ei- 
ther a (sub)giant or a dwarf, depending on the location on 
the original cluster color-magnitude diagram, and the num- 
ber of giants included in the sample was counted. From the 
above calculation we found that the contamination rates are 
at about the 10%— 15% level from the M13 photometry, and 
the 15%-20% level when using the M92 photometry. M92 
is more metal-poor than M13, so more cluster giants fall in 
the color range (0.2 < g — r < 0.5), which is implicitly set by 
our sample selection criteria. This results in a slightly higher 
contamination rate from the M92 photometry. 

More than any other parameters (e.g., shape parameters of 
the halo) in the model, we found that the assumed fraction of 
giants is the dominating factor that determines the total con- 
tamination rate in our sample. Note that our estima te for giant 
contam ination in the sample is highe r than what [ juric et aD 
(2008) estimated (~ 15% vs. ~ 4%). lJuric et all 020081) ob- 
tained an estimate of the fraction of giants to be about 5% 
from M13, using the SDSS pipeline Photo values in this 
crowded region, over a color range similar to that in our anal- 
ysis, and concluded that the bias in the number density of the 
halo is about 4% arising from the misidentification of giants 
as main-sequence dwarfs. 

Giant contamination in our sample produces an overall shift 
of the photometric metallicity estimates toward higher val- 
ues. Figure [14] shows the fiducial sequences on the u - g 
versus g — r diagram for a number of clusters that were used 
in our color calibration. The black lines shown are fiducial 
sequences, with the dotted lines indicating a giant sequence. 
For M67, only a main sequence was detected in SDSS (gi- 
ants are too bright to be included: lAn et alj |2008). The over- 
laid gray lines in each panel show our calibrated models at 
[Fe/H]=-3.0 (left most), -2.4, -1.6, -0.8, and 0.0 (right 
most), respectively. At a given g — r, giants have redder 
u — g colors than the main sequence, which leads to an over- 
estimated photometric metallicity if they are misidentified as 
dwarfs. The size of a bias in the photometric metallicity can 
be as large as 0.5-1 dex, but fortunately this only applies for 
a limited number of stars along each line of sight. 

On the other hand, contamination by thick-disk stars in our 
sample is negligible. We performed a set of Galactic simu- 
lations to check the overall fraction of thick-disk interlopers 
in our sample along various lines of sight, using artificial star 
test results in § 12.51 To simulate a dispersion in the underlying 
[Fe/H] distribution, we combined artificial stars at the central 
[Fe/H] (w -0.7 for the thick disk and « -1.6 for the halo) with 
those at ±ler values (±0.2 and ±0.4 dex for the thick disk and 
the halo, respectively), taking normalizations from a Gaussian 
distribution. We adopted the Galactic structural parameters in 
lJuric et all ( 120081) . 

The fraction of thick-disk stars, of course, is varying at 
different Galactic latitudes, so we computed the fraction of 
thick-disk stars along Stripe 82. We found that the average 
fraction is negligible (0.4% of t he entire sample) b elow pho- 
tometric [Fe/H]= -1.0 (see also IBovv et~aT]|2012bl) . This is 
because our sample selection, based on mass, metallicity, and 
distance, is strongly biased against stars with [Fe/H]> —1.2. 
If all the stars below solar metallicity are included in this esti- 
mate, the fraction becomes ~ 3.4%. However, only stars with 
photometric metallicities less than [Fe/H]=-1 are concerned 
in the following discussions, and we can safely assume that 
the thick-disk contamination is negligible. 

4. RESULTS 
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FIG. 14. — Fiducial sequences of the calibration cluster sample (black lines) in the u — g vs. g — r plane. The solid black lines represent a dwarf sequence; 
dotted lines show the giant sequence. Overlaid gray lines are calibrated models at [Fe/H]= -3.0 (left), -2.4, —1.6, —0.8, and 0.0 (right), respectively. At a given 
g — r, giants have redder u — g colors than the main sequence, leading to an overestimated photometric metallicity. Values in parenthesis next to the cluster name 
indicate the cluster metal abundance in [Fe/H]. 



4. 1 . The Observed Metallicity Distribution Function of the 
Galactic Halo 

Figure [15] shows our in situ observed MDFs from the cali- 
bration (top panel) and the coadded catalogs (bottom panel), 
respectively, includi ng only star s satisf ying the selection crite- 
ria described above. Juri c et all (12008b found that both thick- 
disk and halo stars have approximately the same number den- 
sity at 3 kpc above the Galactic plane. In order to minimize 
the contribution from thick-disk interlopers in our halo sam- 
ple, we imposed a heliocentric distance limit of greater than 
5 kpc and Galactic latitudes \b\ > 35°, which correspond to 
a minimum vertical distance |Z| = 2.9 kpc from the Galac- 
tic plane at the low-latitude limit. As noted above, we im- 
posed a sample cut at [Fe/H] trae = —1.2, where [Fe/JTJtme in- 
dicates a true metallicity value. The great majority of thick- 
disk stars having metallicities above this limit would have 
been excluded from our photometric sample selection. Nev- 
ertheless, it remains possible that some metal-wea k thick- 
disk stars dChiba & Beersll2000t ICarollo et aLlHoTol and ref- 
erences therein) have entered our sample, at least near the 
low-latitude limit. Note that some stars near [Fe/H] true = -1 .2 
could have photometric metallicities greater than [Fe/H] p h t = 
-1 .2, because of non-zero photometric errors. 

The thick solid histograms in Figure Q3] show our derived 
photometric MDFs, cons truct ed using the YREC -based ex- 
tinction coefficients in § 12.31 the dotted line s represent the 
case w ith 10% smaller extinction values than in Schlegel et al. 
(1998). The thin solid histograms are those based o n ex- 
tinction coefficients from Schlafly & Finkbe inerl d201 lb . Al- 
though their extinction coefficients are ~ 20% smaller than 
our default YREC -based values, the difference has only a 



Stripe 82: Calibration Catalog 



Standard 




Stripe 82: Coadded Catalog 




[Fe/H] phot 

FIG. 15. — Photometric MDFs of the halo. The black solid histogram 
shows a MDF using our adopted extinctions and extinction coefficients. The 
dotted histogram is a MDF assuming 10% smaller extinction values; the 
thin solid histogram represen ts a case assuming extinction coefficients from 
Schlaflv & Finkbeiner (201 1). The top panel shows the MDF from the cali- 
bration catalog and the bottom panel shows the MDF from the coadded cata- 
log. The photometric MDFs are complete for [Fe/H] lrl i e < — 1.2. 



Metallicity Distribution Function of the Galactic Halo 



15 



Stripe 82: Calibration Catalog 

: — I 1 1 1 1 1 1 1 1 1 1 1 — : 

- Photometric Metallicity 




Stripe 82: Coadded Catalog 

: — I 1 1 1 1 1 1 1 1 1 1 1 — : 




[FVH] phol 

FIG. 16. — Photometric MDFs of the halo (black solid histograms). The 
top panel shows the MDF from the calibration catalog and the bottom panel 
shows the MDF from the coadded catalog. Our MDFs are complete for 
[Fe/H]nuc < —1.2. The overlaid gray shaded histog rams are MDF s for 
the kinematically-selected local halo stars from Rvan & Norris 1 19 91bfl . for 
which the sample size has been multiplied by a factor of eight, in order to fa- 
cilitate comparison with the photome tric MDFs. The dotted line is the MDF 
from the same set of stars I Ryan & Norris 19 91b ), but after convolving it with 
the photometric convolution kernels (Figure [Tot , in order to simulate the ef- 
fects of photometric errors and unresolved binaries and/or blends (see text). 



moderate effect on the resulting MDFs; the overall shape of 
the MDFs does not change significantly, because of the small 
level of foreground dust toward these high Galactic latitude 
stars. 

The photometric MDF s ar e insensitive to our adopted 
[a/Fe] and age scheme (§ 12.11 ). Even if we had taken a dif- 
ferent [a/Fe] or age relation as a function of [Fe/H], model 
colors would still be forced to match the sa me s et of cluster 
data from our color calibration procedure (§ 12.21 ). The result- 
ing color- [Fe/H] relation would then be essentially unaffected 
by this change. Although intrinsically broad distributions in 
[a/Fe] and/or age at each [Fe/H] could result in additional 
spreading of the observed MDFs, we neglect their effects in 
the following discussion. 

Figure [16] compares our phot ometric MDFs with the spec- 
troscopic MDF obtained by iRyan & Norris! (|T991b, gray 
shaded), which is based on a kinematically-selected sam- 
ple of 372 local halo subdwarfs (c/heiio ^ 400 pc). We have 
multiplied their sample by a factor of eight to directly com- 
pare the shapes of M DFs. The mean [Fe/H] value in the 
IRyan & Norrisl d 199 lb!) sample is -1.80 over -3 < [Fe/H] < 
-1, with a standard deviation of 0.46 dex. Note that, ac- 
cording to these authors, their sample is incomplete above 
[Fe/H]=-1. 

The abundance measurements in IRyan & Norris! <[l991b) 
are accurate to a ~ 0.2 dex, while our photometric metal- 
licity estimates have larger dispersions and asymmetric er- 
ror distributions. In order to perform an "apples-to-apples" 



compari son with our photome tric MDFs, we convolved the 
MDF of Ryan & Norris! dl991bl) with the convolution kernels 
in Figure [10] We binned their data at [Fe/H]< -1 with in- 
tervals of A[Fe/H] = 0.1 dex, and applied the convolution 
kernels to each bin, ignoring spectroscopic [Fe/H] errors in 
IRyan & Norrisl The resulting MDF is shown (with an arbi- 
trary normalization) as a dotted line in Figure [16] which in- 
cludes the effects of both photometric errors and unresolved 
binaries and/or blendin g, as in our photometric MDFs. The 
IRyan & Norrisl (Il991bl) MDFs are similar in their overall 
shape, both before and after the convolution is applied, and 
both match our in situ MDFs rather well. 

As shown in Figure [TO] there is a systematic shift in the 
metallicity scale in the coadded catalog. A Kolmogorov- 
Smirnov (K-S) t est reje cts the null hypothesis that the MDF of 
IRyan & Norrisl (Il991bl) and that of the coadded catalog {bot- 
tom panel) are drawn from the same parent population at sig- 
nificant levels ( p < 0.001). O n the other hand, the null hy- 
pothesis that the lRvan & Norrisl MDF and that of the calibra- 
tion catalog (top panel) are drawn from the same parent pop- 
ulation cannot be rejected (p = 0.364); in other words, their 
MDFs are statistically similar with each other. Photometric 
zero-point shifts may be responsible for the systematically 
hi gher photometric [Fe/ H] esti mates in the coadded catalog 
(§Oand see below in §g221l. 

Perhaps the overall good agreem ent between ou r photo- 
metric MDFs and the spectroscopic IRyan & Norrisl l[l991b) 
MDFs should not be so surprising, given the similar range 
of Galactocentric dist ances of these two samples. The 
IRyan & Norrisl (Il991bl) stars occupy a local volume (although 
their orbits extend to much larger distances), at roughly 
8 kpc away from the Galactic center. Our sample stars in 
Stripe 82 are found at present distances within 7-14 kpc 
from the Galactic center, as shown in Figure [13] although 
they are much farther away from the Galactic disk than the 
IRyan & Norrisl (Il991bl) sample stars. It is worth keeping in 
mind that, although we do not presently have full space mo- 
tions for the stars in our photometric MDFs, they surely con- 
tain many stars with orbits that take them far outside the lo- 
cal volume. Both samples of stars can be thought of as local 
"snapshots" of a volume in the relatively nearby halo. 

4.2. Deconvolution of the Halo Metallicity Distribution 

We now describe our attempts to recover the "true" halo 
MDF, using experiments that consider different ways for un- 
derstanding the nature of the underlying [Fe/H] distributions. 
This requires confronting the possibility that more than one 
parent stellar population may be required, as discussed below. 

We began by deconvolving the observed MDFs of the halo 
using a single [Fe/H] p h t profile (one of the kernels in Fig- 
ure [10), where a single [Fe/H] tlue is assumed, without an in- 
trinsic dispersion. We found that none of these profiles are 
able to match the overall shape of the MDFs, immediately 
suggesting that either the hypothesis of a small dispersion 
in [Fe/H] true should be abandoned, and/or at least two sub- 
components are required to properly account for the overall 
shape of the MDF. We consider each of these two cases in 
turn below. 

4.2.1. A One-Component Model 

Figure [17] shows the fitting results for the observed photo- 
metric MDFs, assuming that the metallicity of the underlying 
halo population can be well-described by a single Gaussian 



16 



An et al. 



Stripe 82: Calibration Catalog 

-r 




Stripe 82: Coadded Catalog 




[Fe/H] 



phot 

FIG. 17. — Deconvolution of photometric MDFs of the halo using a single 
Gaussian [Fe/H] distribution. Histograms are observed MDFs from the cali- 
bration (top) and the coadded catalog (bottom), respectively. The error bars 
represent ±1ct Poisson errors. The solid red line shows an error -con volved 
[Fe/H] p i,ot distribution from the deconvolution kernels in Figure [Tol with a 
peak of the underlying [Fe/H] distribution at [Fe/H]= -1.80 for the calibra- 
tion catalog (top), and a peak at [Fe/H]= —1.55 for the coadded catalog (bot- 
tom), respectively, both with dispersions of 0.4 dex. The dashed blue line is 
the best fitting simple chemical evolution model from Hartwick 1 1976), after 
applying the deconvolution kernels. 



distribution with a dispersion in metallicity. The solid red line 
in Figure [17] shows the resulting fit, obtai ned using the non- 
linear least squares fitting routine MPFIT ( Markwardti l2009l) 
over -2.8 < [Fe/H] p h t < -1.0, after application of the convo- 
lution kernels (FigurefTUb. As described earlier, each of these 
simulated profiles includes the effects of photometric errors 
(oVri = 0.02 mag, <j u z = 0.03 mag) and a 50% unresolved bi- 
nary fraction and/or blends with the M35 mass function for 
secondaries. Note that these kernels are not Gaussian func- 
tions, nor are they the products of a chemical evolution model. 

For the calibration catalog (top panel), we found that a best- 
fitting Gaussian has a peak at [Fe/H] tme = -1.80 ±0.02 with 
® [Fe/H] = 0.41 ±0.03 dex. The reduced \ 2 value (xl) of the 
fit is 2.0 for 15 degrees of freedom, assuming Poisson errors. 
The errors in the above parameters were scaled based on the 
X\, of the fit. The estimated pea k [Fe/H] and the dispers ion are 
similar to those obtained for the lRyan & Norrisl (11991b) sam- 
ple, as expected from the similar MDF shape from these two 
studies. For the coadded catalog (bottom panel), we found a 
best-fitting Gaussian with a peak at [Fe/H] true = -1.55 ±0.03 
and <7 [Fe / H ] = 0.43 ± 0.04 dex, with xl = 2-7. 

For an additional test, we assumed that the photometric 
MDF is shaped primarily by large photometric errors, even 
though the underlying [Fe/H] distribution is single-peaked at 
around [Fe/H]= —1.6. We found marginal agreement with the 
observed MDFs only if the size of photometric errors were 
underestimated by a factor of two (i.e., cr g . r .i w 0.04 mag, 
fj„. z « 0.06 mag) for both of the Stripe 82 catalogs. However, 



we consider it unlikely that the errors have been underesti- 
mated by this much. This hypothesis is also inconsistent with 
the intrinsic ally wide range of spec troscopic [Fe/H] determi- 
nations from lRvan & Norris ( 199 llj. 

The dashed blue line in Figure [T7J show s the best-fitting 
simple chemical evolution model (Hartwick 1976j). We use d 
the mass-loss modified version as in lRvan & Norrisl (|T991b), 
who found an excellent match of this model to their MDF. 
The model MDF is characterized by a single parameter, the 
effective yield (y e fr), which relates to the mass of ejected met- 
als relative to the mass locked in stars. The simple mass-loss 
modified model adopts instantaneous recycling and mixing of 
metal products in a leaky box, and further assumes a zero ini- 
tial metallicity, constant initial mass function, and a fixed ef- 
fective yield. We utilized [Fe/H] in this model as a surrogate 
for the metallicity. After convolving the model MDF with 
the deconvolution kernels (Figure [TOb. we found log 10 y e ff = 
-1.65 ±0.02 (xl = 2.6 for 15 degrees of freedom) for the cal- 
ibration catalog, and log 10 y e ff = -1.37 ± 0.04 (\l = 5.0) for 
the coadded catalog, which simply correspond to the peak of 
the MDF. The goodness of the fit to the calibration catalog 
data is comparable to that using the single G aussian fit, and 
the res ulting effective yield is close to what iRyan & Norrisl 
(Il991bl) obtained from their MDF (logoff = -1 .6). 

4.2.2. A Two-Component Model 

Although a single-peak Gaussian [Fe/H] distribution de- 
scribes the observed photometric MDF rather well, it is not 
a unique solution. Here we consider a two-peak Gaussian 
[Fe/H] distribution fit to the Stripe 82 MDFs. 

The red and blue curves in the top panel of Figure [18] 
show the best matching pair of simulated profiles to the cal- 
ibration catalog MDF, searched using the MPFIT routine 
over -2.8 < [Fe/H] p h t < -1.0. In this fitting exercise, we 
fixed the dispersion of the underlying [Fe/H] distributions to 
® [Fe/H] = 0.30 dex for both Gaussians. The green curve is the 
sum of these individual components, which exhibits an ex- 
cellent fit to the observed profile (xl = 1.9 for 14 degrees of 
freedom). The underlying true [Fe/H] distribution for each 
of these curves is shown in the bottom panel, which exhibits 
peaks at [Fe/H] tlU e = -1.67 ± 0.08 and -2.33 ± 0.30, respec- 
tively. 

The bott om panel in Fig urelTSIalso shows the spectroscopic 
MDF from lRvan & Norrisl (E991b, gray histogram). Our de- 
convolved [Fe/H] distribution matches their observed MDF 
well on the metal-poor side. Their sample is contaminated 
by disk stars above [Fe/H]« -1, while our photometry sam- 
pl e sel ection excludes metal-rich stars with [Fe/H] tl - ue > -1.2 

(§EB. 

The fractional contribution of the low-metallicity compo- 
nent with a peak at [Fe/H] true = -2.33 (the area under the 
blue curve in the top panel of Figure |T8T > to the entire halo 
sample (the area under the green curve in the top panel) is 
24%, with a clearly strong dependence on metallicity. Be- 
low [Fe/H]= -2.0, the contribution from the low-metallicity 
component is 52% of the total numbers of halo stars in this 
metallicity regime. Above [Fe/H]=-1.5, the high-metallicity 
component (the red curve with a peak at [Fe/H] tlue = -1.67) 
contributes 96% of the total numbers of halo stars. 

The relative fraction of the low-metallicity component de- 
pends on our adopted value for the dispersion of the under- 
lying [Fe/H] distribution. In the above exercise, we assumed 
^[Fe/H] = 0.30 dex for each component. However, the contri- 
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FIG. 18. — Top: Deconvolution of the photometric MDF of the halo from 
the calibration catalog, usin g a double Gaussian [Fe/H] distribution. Con- 
volution kernels (Figure [Tol are applied to the underlying [Fe/H] distribu- 
tion to simulate effects of photometric errors and unresolved binaries and/or 
blends in the sample. The red and blue curves show a convolved [Fe/H] p hot 
distribution from each component of the double Gaussian distribution. The 
green curve is the sum of these contributions. The standard deviation of the 
Gaussians are set to <7[F e /H] = 0.30 dex. Bottom: The true underlying [Fe/H] 
distribution, as derived from the [Fe/H] p hot deconvolution shown in the top 
panel. The red and blue curves are the underlying [Fe/H]t rl i e distributions, 
with peaks at [Fe/H]= —1.67 and —2.33, and the green curve is the sum of 
these two. The gray hi stogram shows the observed [Fe/H] distribution from 
Ryan & Norris 1 1991b), normalized to the peak of the distribution. 

bution from the low [Fe/H] component becomes as high as 
42% if the standard deviation is set to cr [Fe/H] = 0-20 dex. On 
the other hand, fixing the value to & [Fe/H] = 0.4 dex leads to a 
similar result as in the previous section, which requires only 
one component in the fit. If o [Fe/H] = 0.2 dex is a reason- 
able lower limit to the dispersion of each of the underlying 
[Fe/H] components, we can only set an upper limit on the 
fractional contribution of the low-metallicity component, at 
~ 40%-50%. 

To investigate this further, we have also performed two- 
component deconvolutions leaving each component's disper- 
sion as a free parameter. F or this, we use the extreme decon- 
volution dBovv et all 1201 11) technique, which fits all param- 
eters of arbitrary numbers of Gaussian components to data 
with individual uncertainties, and which has been extended 
to handle non-Gaussian uncertainty kernels, such as those in 
§ 12.51 The best two-component fit (with Ax 2 = 6.5 compared 
to the fixed o [Fe/H] = 0.3 dex fit above) consists of a 69% 
contribution from a low-metallicity component ([Fe/H]= -2) 
with o [Fe/H] = 0.35 dex and a higher-metallicity component 
( [Fe/H] = -1.4). However, the latter has an unrealistically 
small dispersion of only o [Fe/H] = 0.07 dex. Fitting instead 
with a weak prior on the dispersion (a chi-squared distribu- 
tion with 2 degrees of freedom and a mean of 0.3 dex) yields 
a slightly worse fit (A% 2 = 2 with respect to the best fit), with 
almost equal contributions from a low-metallicity (53%, mean 
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FIG. 19. — Same as in Figure fTsl but for the coadded catalog. The red and 
blue curves in the bottom panel show the underlying [Fe/H] lrlle distributions 
with peaks at [Fe/H]= —1.39 and —1.99, respectively. 

[Fe/H]= -2.14, a [Fe/H] =0.3 dex) and a higher-metallicity 
(mean [Fe/H]=-1.5, u [Fe/H] = 0.2 dex) component. These re- 
sults further show that the parameters of a two-component fit 
are only poorly constrained on the basis of the photometric - 
metallicity data alone. 

An even stronger constraint on the fraction of the low- 
metallicity halo component can be obtained if our de- 
convolved MDF is f orced to match the shape of the 
iRvan & Norris! d!991bl) MDF. As shown in the bottom panel 
of Figure QjH our assumed [Fe/H] dispersion of 0.3 dex re- 
sults in a good agreement of our MDF with that of the Ryan 
& Norris sample (gray shaded). However, the agreement 
breaks down when a smaller (or larger) value is assumed 
for a dispersion of the individual underlying [Fe/H] distri- 
butions; a smaller dispersion leads to an enhanced contribu- 
tion from the metal-poor component, resulting in an overes- 
timated number of metal-poor stars with respect to the Ryan 
& Norris MDF. From the above considerations, we find that 
individual components of the underlying MDF with ~ 0.25- 
0.32 dex dispersion in [ Fe/Hl provide a reasonably good 
fit to the IRvan & NorrTsI (Il991bl) MDF; in other words, the 
Rvan & Norris MDF is also well fit by double Gaussians with 
these dispersion values. The fractions of the low-metallicity 
component of the halo corresponding to these values are 
~ 20%-35% of the entire halo population. 

In the case of the coadded catalog (Figure [19), the over- 
all metallicities are shifted toward higher [Fe/H] values; peak 
[Fe/H] tlue values of the two components are [Fe/H] = -1 .39 ± 
0.12 and -1.99 ±0.23, respectively, with a total xl value 
of 2.6. The shift in [Fe/H] is expected from the photo- 
metric z ero- point difference between the two Stripe 82 cat- 
alogs (§ 13.11 ), where the larger (fainter) M-band magnitudes in 
the coadded catalog produce higher photometric metallicities 
than in the calibration catalog. The contribution of the metal- 
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poor component is comparable (30%) to that derived from the 
calibration catalog (24%). 

A strong disagreement with the iRvan & Norrisl ([l991b) 
MDF is evident in the lower panel, where the deconvolved 
MDF shows systematically higher metallicities. If we expect 
both MDFs to be the same at similar Galactocentric distances, 
this would indicate that the MDF in Figure[T9]has been shifted 
toward higher [Fe/H] d ue to photometric zero-point errors in 
the coadded catalog (§ 13. U . Some of the stars in Stripe 82 
have spectroscopic estimates from the SSPP, and were in- 
cluded in the comparison of the spect rosc opic metallicities 
with the photometric determinations (§ 12.41 ). However, there 
were not a s uffic ient number of stars sati sfying both the spec- 
troscopic (§ 12.41 ) and photometric (§ 13.11 ) selection criteria for 
a meaningful comparison to verify which photometric catalog 
provides a more consistent result with the SSPP estimates. 

Our results illustrate a simple limiting case, where one 
could imagine that the stellar halo formed out of individual 
subhalos which, in the mean, were pre-enriched in metals that 
peak roughly at the observed locations. Nevertheless, the fact 
that the observed photometric [Fe/H] distribution is well de- 
scribed by primarily two metallicity peaks at [Fe/H|~ -1-7 
and —2.3 is consistent with the argument by Carollo et al. 
(120071 120101) that our Milky Way stellar halo is a superpo- 
sition of at least two components, the inner and outer halos, 
that are distinct in metallicity, kinematics, and spatial distri- 
butions. The apparent duality of the Galactic halo system is 
further supported by a limited kinematic anal ysis of stars in 
our photometric sample, as described below (§ 14. 3b . It should 
be kept in mind that our description of these two populations, 
following Carollo et al., assumes that they are "smooth" dis- 
tributions, even though it is recognized that there exist numer- 
ous, possibly even locally dominant, examples of substructure 
(such as stellar streams and overdensities) present as well. 

From the above tests, we conclude that we do not have 
sufficiently strong constraints, from the observed photometric 
MDFs alone, in order to discriminate between the hypothesis 
of the duality of the halo (Figures [T8~l - [T9T l and the single stel- 
lar population (Figure [TVj model. We also cannot preclude 
the possibility that the halo comprises more than two sub- 
components, because of the limited number of constraints pre- 
sented by the observed MDFs to perform multi-component fit- 
ting. However, we show below that the degeneracy in this so- 
lution can be at least partially removed by combining indepen- 
dent information from the avail able kinematics of our s ample 
stars, as well as for those in the IRvan & Norrisl (11991b!) sam- 
ple. 



4.3. Metallicity Correlations with Kinematics 

Most of the stars in our Stripe 82 photometric sample do 
not have full velocity information available. However, at 
high Galactic latitudes, proper-motion measurements alone 
can provide useful constraints on the U (towards the Galactic 
center) and V (in the direction of Galactic disk rotation) veloc- 
ity components, which are sufficient to constrain the motions 
of stars parallel to the Galactic plane. 

We searched fo r proper-motion measurements in SDSS 
dMunn et al.1 120041) . with the proper-motion flags set to the 
recommended values: match = 1, nFit = 6 and dist22 > 7, 
sigRa < 525, and sigDec < 525 (taking into account the cor- 
rection described by Munn et al. 2008). We then combined 
proper-motion measurements with our distance estimates to 
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FIG. 20. — Distribution of approximate rotation velocity in the Galactocen- 
tric cylindrical system (va) for the calibration (left panels) and the coadded 
(right panels) catalogs. The thick black and thin red histograms are those 
based on stars at \b\ > 45° and > 60°, respectively. The top to bottom 
panels show the sample in different metallicity bins. A vertical dotted line is 
shown at zero rotation. 

compute approximate rotation v elocities (v^) in the Galac- 
tocentric cylindrical system (see Bond et al.H2Q10t) . 

Figure|20]shows the approximate velocity distribution in 
for the calibration and coadded catalogs, respectively. The 
top panels show the full velocity distribution of all stars in 
each catalog, while the middle and bottom panels show sub- 
sets of these in the metal -rich (—1.5 < [Fe/H] < —1.0; middle 
panel) and in the metal-poor (-3.0 < [Fe/H] < -2.0; bottom 
panel) regimes, respectively. Inspection of the distributions in 
these two metallicity ranges indicates that they are not iden- 
tical, as might be expected for a single halo population. The 
thick black histograms show the distributions of stars at 
\b\ > 45°; red histograms show the more restricted case with 
\b\ > 60° (the restriction to higher latitude provides a better 
approximation for the rotational velocity). Comparisons be- 
tween the black and red histograms show that the latitudinal 
restrictions only mildly change the overall distributions in all 
metallicity bins. For the calibration catalog, a K-S test re- 
jects the null hypothesis that the samples in the two different 
metallicity ranges are drawn from the same parent population 
at significant levels (p = 0.019 for \b\ > 45°; p = 0.003 for 
\b\ > 60°). For the coadded catalog, the K-S test is unable to 
reject this hypothesis (p = 0.167 for \b\ > 45°; p = 0.297 for 
\b\ > 60°). Although not a strong effect, there exist greater 
fractions of more metal-poor halo stars in retrograde, rather 
than prograde, rotation for both catalogs. 

A more striking difference is obtained when the MDFs are 
considered for different ranges, as shown in Figure [2T| As 
before, the black and red histograms indicate the two sam- 
ples at different Galactic latitude ranges. The middle pan- 
els show the MDFs for stars with high prograde rotation 
(v^ > 80 km/s); bottom panels show those of stars with ex- 
treme retrograde rotation (v^ < -80 km/s). Note that the 
±80 km/s limit was set in order to maximize the numbers 
of stars in each split, so as to enable statistically meaningful 

16 Below we use to represent a rotational velocity computed using the 
full kinematic information, to distinguish it from va . 
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FIG. 21 . — MDFs in different bins of approximate rotational velocity in the 
Galactocentric cylindrical system (v<^) for the calibration (left panels) and 
the coadded (right panels) catalogs. The thick black and thin red histograms 
are those based on stars at \b\ > 45° and \b\ > 60°, respectively. A vertical 
dotted line at [Fe/H] = -1.6 is shown for reference. 



interpretations. 

Inspection of Figure[2T]indicates that the MDFs for stars in 
retrograde motion are shifted toward lower metallicities than 
the MDFs for those in prograde rotation. For the calibra- 
tion catalog, whether the samples are restricted to \b\ > 45° 
or \b\ > 60 , a K-S test rejects the null hypothesis that the 
two samples for different ranges of are drawn from the 
same parent population at highly significant levels (p = 0.0001 
for \b\ > 45°; p = 0.0002 for \b\ > 60°). For the coadded 
catalog, the lower-latitude cut results in a significant rejec- 
tion (p = 0.026 for \b\ > 45°); no rejection can be made for 
\b\ > 60° (p = 0.257). We conclude that greater fractions of 
low-metallicity halo stars exist in retrograde, rather than pro- 
grade, rotation. 

Inspired by these results, we decided to check how the 
MDFs of the spectroscopic iRyan & Norrisl (Il991bl) sample 
vary for prograde vs. retrograde rotation, and how the 
distributions vary for higher- and lower-metallicity divisions. 
Note that we distinguish from v$ to indicate the fact that 
the rotation component for the Ryan & Norris sample is ob- 
tained from full space motions (iRvan & Norris|[1991al) . Fig- 
ure[22]shows the [Fe/H] distributions (left panels) and dis- 
tributions (right panels) of their local kinematically-selected 
sample, with the same limits as on and [Fe/H] in Fig- 
ures l20l2Tl except in the middle left panel, where we set 
Vcf, > 50 km/s to include a sufficient number of stars. There 
exist clear differences in these distributions. A K-S test in- 
dicates that the null hypothesis that the samples are drawn 
from the same parent population is rejected at significant lev- 
els: p = 0.019 for the different ranges in and p = 0.021 
for the different ranges in [Fe/H]. Both divisions indicate that 
greater fractions of metal-poor halo stars are found in retro- 
grade, rather than prograde, rotation. 

Table [3] summarizes the average [Fe/H] values (and errors 
in the mean) for the kinematically-divided samples consid- 
ered above. For the calibration and coadded catalogs, the 
mean [Fe/H] is between —1.7 and —1.5 for stars in highly 
prograde motion (v$ > 80 km/s). For the calibration cata- 
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FIG. 22. — Same as in Figures 120121 1 but for the local kinematically- 
selected sample of Rvan & Norris 1 1991b). Note that we distinguish from 
to indicate the fact that the rotation component for the Ryan & Norris sam- 
ple is obtained from full space motions. The vertical dotted line is shown in 
the left panels for [Fe/H] = —1.6, and at zero rotation in the right panels, for 
reference. 

log, the mean [Fe/H] for the highly retrograde components 
(v^ < -80 km/s) is ~ 0.3-0.4 dex lower than for the metal- 
rich counterparts; the difference is ~ 0.2 dex for the coad- 
ded catalog. The difference in the mean [Fe/H] between the 
highly retr ograde component and th e highly prograde compo- 
nent in the Rvan & Norris ( 1991b) sample is 0.3 dex. Taken 
as a whole, these simple kinematic divisions indicate that a 
single stellar population is inadequate to describe the behav- 
ior of the higher and lower-metallicity halo stars for all three 
samples. 

We caution that the observed kinematic-metallicity correla- 
tion can also be produce d by (sub)giant contamination in our 
sample (e.g., Ryan 1989). Unrecognized giants and subgiants 
are brighter than dwarfs, so their distances and velocities will 
be underestimated, while their photometric metallicities are 
higher than those for dwarfs (F igure [T4l. We used the pho- 
tometry of Ml 3 (lAn et al.l 12008b to per form the same suite 
of simple Galactic halo models as in § 13.21 and found that 
(sub)giants have systematically higher photometric metallic- 
ities, and higher v^, than dwarfs. The overall contamination 
rate of (sub)giants in the sample is estimated around 15%, but 
the fraction increases to ~ 40% when high-metallicity stars 
(-1.5 < [Fe/H] < -1.0) are considered. However, it is less 
likely that the observed kinematic-metallicity correlation is 
caused entirely by the (sub)giant contamination, because the 
IRyan & Norrisl sample, which is independent from our pho- 
tometric sample and has different selection criteria, also ex- 
hibits the same behavior (Figure [22). Further investigation is 
needed on the effect of (sub)giant contamination. 

5. SUMMARY AND DISCUSSION 

We have used a new, empirically calibrated set of stellar 
isochrones to derive distances, temperatures, and metallici- 
ties for individual stars in SDSS Stripe 82 using ugriz pho- 
tometry. These estimates have been validated by comparisons 
with the IRFM temperature scale and SDSS/SEGUE medium- 
resolution spectroscopic values. Based on unbiased photo- 
metric samples of main-sequence Stripe 82 stars in a relatively 
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TABLE 3 

Mean [Fe/H] of Kinematically Separated Samples 



Latitudinal 


> 80 km/s 


V4, < -80 km/s 


Sample 


restriction 


(LFe/HJ) 


N 


(LFe/HJ) 


N 


Stripe 82 Calibration 


b 


> 


60° 


-1.69 ±0.05 


89 


-2.08 ±0.08 


41 


Stripe 82 Calibration 


b 


> 


45° 


-1.72 ±0.04 


143 


-2.05 ±0.05 


75 


Stripe 82 Coadded 


b 


> 


60° 


-1.47 ±0.05 


59 


-1.63 ±0.06 


56 


Stripe 82 Coadded 


b 


> 


45° 


-1.51 ±0.04 


102 


-1.69 ±0.05 


111 


Rvan& Norris (^lbV 1 








-1.59±0.09 b 


47 


-1.90 ±0.07 


88 



b Rotation component for the Ryan & Norris sample obtained from full space motions 
(Vj,) is used to divide the sample in the above table. 
b V<4 > 50 km/s was used. 



narrow mass range, we have constructed an in situ MDF of 
the Galactic halo in the distance range 5-8 kpc from the Sun, 
which is similar in s hape to that of t he loca l kinematically- 
selected subdwarfs in Ryan & Norris (1991b). This MDF can 
be adequately fit using a single Gaussian with peak at [Fe/H] 
= -1 .5 and a broad (0.4 dex) dispersion of metallicities. How- 
ever, an equally valid fit can be obtained from the use of 
two components with peaks at [Fe/H] rs -1.7 and -2.3, with 
smaller dispersions (0.25-0.32 dex). For an adopted disper- 
sion of 0.3 dex, on the order of 25% of the stars in our local 
halo sample can be associated with the low-metallicity com- 
ponent. For stars with [Fe/H] < -2.0, this fraction increases 
to on the order of 50%. For stars with [Fe/H] > —1.5, essen- 
tially all of the stars are associated with the high-metallicity 
component. If a smaller dispersion is adopted (0.2 dex), 
~ 40% - 50% of the total halo sample can be associated with 
the low-metallicity component, but would make the decon- 
volv ed photometric MDF dif ferent in shape from that of the 
local iRyan & Norrisl (Il991bl) kinematically-selected sample. 
Our analysis cannot preclude the possibility that more than 
two sub-components could exist in the halo. 

A limited kinematic analysis of the stars in our photometric 
sample, based on approximation of the rotational motions of 
stars about the Galactic center (derived from available proper 
motions at high Galactic latitude), indicates that the metal- 
poor and metal-rich subsamples exhibit different behaviors, 
suggesting that the dual-component model is likely to be su- 
perior to the single-component model. Greater fractions of 
metal-poor stars ([Fe/H] < -2.0) in the halo possess retro- 
grade orbits, as compared to more metal-rich stars (-1.5 < 
Fe/H] < — 1 .0), consistent with the claim by Caro llo et al.1 
20071 120101) that the halo comprises two overlapping sys- 
tems, an inner-halo population and an outer-halo popula- 
tion, with similar metallicity-kinematic correlations as we 
find from our present analysis. The observed behaviors of 
these halo stars are also consistent with theoretical expec- 
tations from Milky Way-like galaxies, where inner-spheroid 
stars primarily formed in the main progenitor(s) of the galaxy 
and exhibit net prograde rotation, while the accreted outer- 
halo component exhibits a mild degr ee of rotation, or in a m i- 
nority of cases, a retrograde rotation (iMcCarfhy et al.ll2012l) . 

We have also considered an independent sample of halo 
stars, one that is still in co mmon use for describi ng the 
MDF of the Milky Way's halo. IRyan & Norrisl dl991bl) found 
that their derived halo MDF (based on a local kinematically- 
selected sample) can be well-describ ed by a simple chemi- 
cal evolution model (Hartwickl ll976l) . Our local in situ halo 
MDF is also well-fit by this model. However, there are a 
number of unphysical assumptions that are made by such a 



model, including instantaneous recycling and mixing of met- 
als throughout the entire halo, as well as the use of a so-called 
effective yield (by which the metals produced by individual 
supernovae are reduced by an arbitrary amount, one that is not 
predicted by present understanding, but acts only as a fitting 
parameter of the model). Such limitations immediately call 
into question the applicability of these classes of models in 
the context of a presumed hierarchical assembly of the Galac- 
tic halo from individual lower-mass subhaloes. Such a model 
(at least when applied to the halo as a whole) clearly presumes 
that the stars belong to a single population. As was the case 
for our analysis of the unbiased photometric stellar halo sam- 
ple, our investigation of the rotational behavior of the Ryan & 
Norris data also does not support a single-population model. 

Our current understanding on the nature of the halo pro- 
vides caveats for the use and interpret a tion of the ob served 
MDFs . For example, ISchorck et all (120091) and iLi et all 
(2010) consider halo MDFs based on spectroscopic follow- 
up of candidate metal-poor stars from the Hamburg/ESO sur- 
vey. Below [Fe/H]= -2.5, the MDFs of both samples (the 
first comprising mainly giants, the second comprising mostly 
dwarfs) appear in excellent agreement with one another, and 
led to the conclusion that there existed a sharp low-metallicity 
cutoff (at [Fe/H]= -3.6), below which fewer than expected 
numbers of the lowest metallicity stars were found. For 
stars with higher metallicity, these authors argue that the ob- 
served MDFs of the two samples differ from one another due 
to the very different selection biases on metallicity involved 
for the cooler (giants) and warmer (dwarfs) samples. How- 
ever, i f the halo comprises two overlapping stella r compo- 
nents dCarollo etalJ 120071 l2010t iBeers et aljl2012h . the dif- 
ferent behaviors of the giant and dwarf samples are likely to 
be strongly influenced by the degree to which the inner- and 
outer-halo populations contribute to them, convolved with the 
clearly present selection biases, rather than being due to the 
selection bias alone. The implication of the dual nature of 
the halo is that the subject under discussion is not the tail of 
the MDF of "the halo," but rather, the tail of the MDF of the 
outer-halo population, with clear implications for the origin of 
this component. Further investigations of these connections 
would be of great interest. 
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